Wetting transition on a one— dimensional disorder 



D.M. Gangardt 1 , S.K. Nechaev 2 * 
1 School of Physics and Astronomy, University of Birmingham, Edgbaston, Birmingham B15 2TT, UK 
2 LPTMS, Universite Paris Sud, 91405 Orsay Cedex, France 

(Dated: February 5, 2008) 

We consider wetting of a one-dimensional random walk on a half-line x > in a short-ranged 

■ potential located at the origin x — 0. We demonstrate explicitly how the presence of a quenched 
| chemical disorder affects the pinning-depinning transition point. For small disorders we develop a 

perturbative technique which enables us to compute explicitly the averaged temperature (energy) 
£SJ ' of the pinning transition. For strong disorder we compute the transition point both numerically 

and using the renormalization group approach. Our consideration is based on the following idea: 
fl P the random potential can be viewed as a periodic potential with the period n in the limit n — > oo. 

■ The advantage of our approach stems from the ability to integrate exactly over all spatial degrees 
of freedoms in the model and to reduce the initial problem to the analysis of eigenvalues and 

04 , eigenfunctions of some special non-Hermitian random matrix with disorder-dependent diagonal and 

■ constant off-diagonal coefficients. We show that even for strong disorder the shift of the averaged 
pinning point of the random walk in the ensemble of random realizations of substrate disorder is 

'(—{'', indistinguishable from the pinning point of the system with preaveraged (i.e. annealed) Boltzmann 

. weight. 

B 

J_> , I. INTRODUCTION 

a • 
•i— > . 

Wetting is one of the most intensively studied phenomena of statistical physics of interfaces. In a very general 
& ' setting wetting implies the interface pinning by a solid impenetrable substrate. Problems of interface statistics in the 
presence of a hard wall were addressed in many publications (see, for example, [l[ and references therein). The most 
I interesting questing concerns the nature of the wetting or pinning-depinning transition of the interface controlled by 
parameters of its interactions with the substrate. Here we study the case when the substrate is inhomogeneous, so 
^ , the wetting transition occurs in the presence of disorder. 

O . The pinning-depinning transition in models of wetting in presence of quenched disorder was studied by many 
research groups since the middle of 80th. In 1986 Forgacs et al Q developed a perturbative renormalization group 
r- — i ■ approach to the (l+l)-dimensional wetting subject to a disordered potential along the substrate. Around the same 
^ ' time Grosberg and Schakhnovich Q applied the RG technique for studying an equivalent problem of the localization 
C*") , transition in ideal heteropolymer chains with quenched random chemical (primary) structure at a point-like potential 
ON ■ well in a D-dimensional space. Many conclusions of Q for D=3 agree with those of @- Both approaches provide 
| important information about the thermodynamics near the point of transition from delocalized (depinned) to localized 
(pinned) regimes in the presence of quenched chemical disorder. 

f^S ' However some crucial questions of pinning-depinning transition in a quenched random potential still remains open. 
, One of the most intriguing problems is the determination of the averaged transition temperature, T q , for quenched 
' chemical disorder. Since the temperature enters into problem through the Boltzmann weight (3 = e u "'^ T , where u m 
, is the energy of m-th interface segment, one may attempt to relate the transition point to the temperature T a for 

v~j ■ annealed chemical disorder with preaveraged Boltzmann weight, (/?) = (e ll '"/ T ). The RG approaches 0, Q claim 
T q = T a in the thermodynamic limit. In 1992 Derrida, Hakim and Vannimenus Q have reconsidered the (1+1)- 

■ dimensional model of wetting and have shown by a different RG technique that the disorder is marginally relevant, 
i.e. any infinitely small disorder displaces the averaged transition point T q in the ensemble of quenched sequences from 
the transitionpoint T a in ensemble of sequences with initially preaveraged (annealed) Boltzmann weight. Subsequently 
other works [fj, 1|| arrived at the same conclusion. The equivalent problem of localization transition of a random walk 
has been also deeply studied in mathematical literature. In Q it was rigorously proven that in systems with return 
probability which scales as ~ N~ a , where N is the number of steps, the phase transition curves for quenched and 
annealed systems coincide for 1 < a < 3/2 for small disorder and are different for 3/2 < a < 2 though they are very 
close numerically (again for small disorder). The similar conclusion has been drawn in the work Q by an alternative 
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method. In other work [9( it has been shown rigorously that the disorder is marginally relevant for a > 3/2. As for 
the case a = 3/2 there is no definite answer (even for small disorder) whether the results for quenched and annealed 
disorder coincide. The value a = 3/2 of the critical exponent considered here is, therefore, of particular interest. 

In our paper we demonstrate explicitly how the presence of quenched chemical disorder affects the pinning— depinning 
transition point. For small disorder we develop a perturbation theory which enables us to compute explicitly the 
transition temperature of the system. The advantage of our approach, which borrows the basic idea from (Toj |. stems 
from the ability to integrate exactly over all N spatial degrees of freedom of the model and to reduce the initial 
problem to the analysis of eigenvalues and eigenfunctions of some special non-Hermit ian N x N random matrix with 
disorder-dependent diagonal and constant off-diagonal elements. Our approach is based on the following general idea: 
the random potential can be viewed as a periodic potential with the period N in the limit N — > oo. 

The paper is organized as follows. In Section|TT]we describe the model and derive basic equations. In Section HTT1 we 
develop the perturbation theory for eigenvalues of our random matrix and derive the corresponding expressions for the 
averaged temperature of pinning- depinning transition for any type of disorder (not necessary to be Gaussian). The 
simple renormalization group approach is developed in Section IIVI while the numerical analysis of general analytic 
equations is performed in Section [VJ In Conclusion we summarize our results and pose some new questions. 



II. RANDOM MATRIX FORMULATION OF WETTING PROBLEM 



The problem of fluctuating interface in thermodynamic equilibrium maps into an equivalent problem one- 
dimensional random walk on a half-line x m > 0, where m is the discrete time. On the line x = additional 
Boltzmann weights j3 m = e u " ll>T account for the potential interaction of the fluctuating interface (the random walk) 
and the impenetrable substrate. As soon as the running time t = m can be associated with the current coordinate 
along the wall, one can say that u m is the interaction energy of an interface with a substrate at a position m. If the 
interaction energies u m are arbitrary, then the set {u m } = {ui, u%, u^} represents the quenched random interaction 
of the interface with the substrate. 

Consider a random walk in a semi-axis x > which represents the height of the fluctuating interface interacting 
randomly with a surface situated at x = 1. The probability of the random walk interacting with random surface 
potential {u m } to be found at the position x after TV steps will be denoted by Gn(x). This function satisfies the 
following recursion relation 

' G N+1 (x) = ±G N (x-l) + lG N (x + l) + ±(p N -l)5 x , 1 G N (x + l) x>l 
< G N (x)=0 x = (1) 

Gn=o{x) = 5 x ,i 

In the presence of the disordered potential u m ^ the equation (TTJ does not conserve normalization of the propagator 
Gjv(x) so it has to be explicitly normalized after N steps. The typical configuration of the random surface is depicted 
in Fig[Tjfor bimodal disorder {u m } = uo or u\. Let us stress however that our considerations are quite general and are 
not restricted to any specific type of a substrate disorder. To answer the question about the location of the pinning 




FIG. 1: Wetting in a random potential on a line. 

transition find it is more convenient to change distribution of (3 m by changing parameters of u m at fixed temperature 
T which is conventionally set equal to unity for the rest of the article. 
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Now we define explicitly what is the pinning (or localization) of the random interface in (l+l)-dimensional wetting 
problem. Consider the mean-square end-to-end distance, (cc 2 (iV)) of the random interface of length N 



(2) 



There exists some critical value of the energy, u tr = ln/3 tr , which separates two different types of behavior of the 
mean-square end-to-end distance in the thermodynamic limit: 



Ai(u) for u > u tI 
A2(u)N for u < Ut T 



(3) 



where Ai^{u) are some positive energy-dependent constants independent of N. The value Ut T is called the energy 
of the pinning transition. Formally speaking, for all u > u tl - the Eq.JT]) has a discrete spectrum, and for u < u tT 
Eq.(fT]) has a continuous one. Slightly above the transition point the following critical behavior of the free energy 
-F/v (u) = In Gn of an A'-step random walk in a half-space x > is expected: 



f(u) = lim ^F N {u) 

N—>oc 1\ 



lim — In Gn 

N^oc N 



const (u — Uti) a 



(4) 



where a is the critical exponent defining the order of the phase transition. In the absence of disorder it was shown 
that a = 2 which corresponds to an ordinary 2nd order phase transition. In the next subsection we reproduce 
these calculations. The reason for doing this is to define our notation and introduce important concepts which we 
shall use in the more complex disordered case. 



A. Wetting in absence of substrate disorder 

We review briefly the situation with no chemical disorder , i.e. when all interaction energies take the same value, 
{u m } = {u, u, . . . , u}. Hence, all the Boltzmann weights in Eq.{T]) are equal, i.e. (3 m = f3 for m = 1, ...N . 

Using the discrete sin-Fourier transform and introducing the generating function, we define the function G(q, s) as 
follows 



G(q, s) = sin go; s n Gn{x) 



x=0 



N=0 



This function satisfies the following integral equation 



Introducing 



—G(q,s) ■ — - = cos q G(q, s) H — (8— l)sin<7 / G(q' ,s) sin 2q'dq' 

S S 7T J 



(G) 



A(s) = / G(q,s) sin 2qdq 



we can rewrite (|6]) as an algebraic equation for A(s 

r7V sin q sin 2q 



A(s) 



dq + -{B-l)A(s) —2 Uq 

1 — s cos q tt J 1 — s cos q 



Solving it we get 



A(s) = 



1 



7T 



(i - Vi^f 



sin q sin 2q 
o 1 — s cos q 



dq 



s s s 2 

--G9-1) 

TT 



(7) 



(8) 



(9) 
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which allows us to write the complete expression for the function G(q, s): 



1 — s cos q 

The inverse Fourier transform applied to (|10[) gives us 



G( q ,s) = T ^(l + S -^A { s)) (10) 



1 (3- 1 . . A / 1 -y/l 



s 



2 



G(x,s) = - G(q,s)sinqxdq = 2[- + - A(s)){ . (11) 

7T J \S IT /\ S / 

Performing the summation over all x > we arrive at the following expressions for the functions G(s) ~ X)^Lo G(x, s ) 
an d J2T=o x2 G(x, s): 



G(s)=J2G(x, S )=2(i + ^J1a(s: 



£ x*G(x, s) = 2 ( 1 + '^1A( S) ^ ^ + ( 2 + 'J'" 1 + 



a=0 



(12) 



(s - 1 + VT^72)2 



The first equation in (|12[) is the basis for the derivation of the free energy / = lirn/v^oo In Gat and we have 

G N = ^-j>G{s)s- N - 1 ds, (13) 
while the second equation in (|12p provides the mean-square end-to-end distance 

(x 2 (N)) = G N ' ^- (f s-"- 1 dsjTx 2 G(x,s). (14) 

The integrals in Eqs. ()13ll4p are performed along a contour in the complex plane of s which lies inside the unit circle. 
The function G(s) is analytic in the whole complex plane except for the branch cuts on the real axis for |Res| > 1. In 
the localized phase the function G(s) has poles for |s| < 1 on the real axis. In the thermodynamic limit N — ► oo the 
pinning transition is indicated by appearance of the pole at s — > s tr — 1, i.e. is determined by the point of divergence 
of the function G(s) for s = 1, which, in turn, diverges when the denominator of A(l) tends to zero as (3 approaches 
the transition point (3t r = e" tr - Thus, we have the following equation for /3 tr (see Eq.Q) 



Str 



1- a/1-4 



= 0; p tI = e Ut '=2; u tI = hx2. (15) 

=i 



In the localized phase the free energy / in the thermodynamic limit N — > oo is dominated by the closest to zero 
pole sq((3) of the function G(s) in (|12p for some fixed value of (3: 



1=0; so = 2 ^^ - (16) 



Calculating the contribution of this pole to the integral in Eq. (|13[) . so that: 

1 I _-N-l > _ „-JV f 2 Vfl - ! X 



G - = A, 2™ / = °o = [^T- ) (17) 



we arrive at desired expression of the free energy / in the thermodynamic limit 



/ = -ln S0 =ln^O. (18) 

P 

In the vicinity of the transition point we write (3 = f3 tr + S, where S <Si 0tr- Expanding (|18|) near [3 tT = 2, we get 

Comparing (flU)) to ^ we conclude that a = 2 and hence the pinning transition on a homogeneous substrate is the 
standard 2nd order phase transition in accordance with the results 
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B. Wetting in a periodic bimodal potential 



Consider now the wetting in a substrate potential with a bimodal periodic distribution of energies {«} = 
{iio, tii, Uo, iti, uo, ui}. We denote by uq and u\ the energies belonging to the even/odd time slices correspond- 
ingly. This problem has been addressed for the first time in [12j and then considered in much more general setting in 
subsequent publications [HI, 0, EE EH ■ The reason to reconsider this problem is basically methodological: we solve 
this problem in a matrix form and then in Section III CI generalize this matrix approach to a substrate with arbitrary 
period of disorder. 

The master equation for the function Gn(x) which generalizes ([T]) is as follows: 

G 2N+1 (x) = ±G 2N {x - 1) + \G 2N (x + ].) + !(#,- l)5 xA G 2N (x + 1) 
G 2N+2 {x) = \G 2N+ i{x - 1) + \G 2N+ i{x + 1) + \{(3i - l)S X: iG 2N+ i(x + 1) 
G N (x = 0) = 
. G N =o(x) = 5 x> i 

where /3q,i = e" ' 1 are the corresponding Boltzmann weights. 
Define odd and even functions Gn- 



G 2N {x) = W N (x) 
G 2N +i(x) = V N (x) 

Rewrite (|2T))) in Fourier space using functions Wn(x) and Vn(x): 

V N {q) = COSqW N (q) + — (/3 - 1) / W N (q') sin 2q' dq' 

I" JO 

sin a 2 f" 

W N+1 {q) = COS qV N (q) + —± / V N (q') sin 2q' dq' 

n n Jo 

, W N=0 (q) = sinq 
Introducing the generating functions W(q, s) and V(q, s) 



(20) 



(21) 



(22) 



W(q, s) = £ W(q)s N ; V(q, s) = £ V{q)i 



N 



N=0 



N=0 



we can write a closed system of integral equations 



V(q, s) = cos q W(q, s) + 



sing / W(q', s) sm2q'dq' 



W(q, s) = sin q + s cos q V(q, s) + s 



j3\ 1 /" 7r 

sing / V(q', s) sin 2q'dq' 

7T 7 



(23) 



(24) 



It is worth noting that the variable s plays the role of the fugacity of the two-step block and is no longer associated 
with a single step as, for example, in Eq.([T]). 

Equations (|24[) allow for a very convenient matrix formulation which could be later easily generalized to longer 
periods. Introduce the matrices 



(25) 











: 











1 
















c 


: 













and vectors 



G(q,s) 



V(q, S) 
W(q,s) 



F(g) 




sing 



(26) 
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Rewriting Eq. ([24|) using (j25]) -(|27 [) . we obtain equation for the vector function G(q, s) 



sing 



G(q,s) = F(g) + cosgM s G(g,s) 



For further analysis it is convenient to rewrite (|27p in the following form 



BM a / G (q',s) sin 2q' dq' 



Define now 



G(g,s) = [I-cosqMA F(q) 



sing 



J-cosgM s ) BM S \ G(q' , s) sin 2q' dq 



G(g, s) sin 2q dq 



(compare to Eq.([7])). The solution for Q(s) reads 

Q00 



[ dq sin q sin 2q (i — cos q M s \ BM S 



I — cosqM s ) F(g)sin2g<ig 



(27) 



(28) 



(29) 



(30) 



This equation extends the solution ([9]) to the periodic bimodal potential {u} = {uq,u\,uq. > Vj\,...,Uo,Ui\. Setting 
parameter s to its critical value s = Str = 1 and calculating explicitly the integrals we arrive at the following result 



Q(«) = ( V 1 



o 

i 

2-/3o 



2- 




(31) 



implying the phase transition at [3\ = 2. It is clear from the above expression that the value of /3q is actually irrelevant. 
This fact reflects the peculiarity of the microscopic model: after an even number of steps the random walk has exactly 
zero probability to reach x = 2 from which transition to x = 1 is controlled by (3q. Mathematically it is reflected in 
the orthogonality of F(q) to the eigenvector belonging to the eigenvalue 1/(2 — (3 ) in Eq. (|3"Tj) . We shall see in the 
next Subsection that this peculiarity persists for arbitrary even periods of the substrate potential. 



C. Wetting in a potential with arbitrary period length 



We can straightforwardly generalize the approach developed in the previous Section to the case of a substrate 
potential with the period n: 



period period 
{u} = {. . .,U ,Ui, ...,U n -i, ...,U ,Ui, ...,U n -i, . . .} 



(32) 



i.e. the total substrate consists of £ = N/n copies of random subchains of length n each. The equations (written 
already in the Fourier space) which extend (f2"2"| to the case of repeating n-periodic potential are as follows 

G${q) = cos qG%\q) + ^(,3 Q -l) f (q') sin 2q' dq' 



G%\q) = cos qG^(q) + ^% - 1) f G$ {q') sin 2q> dq' 

n Ja 



G%- 1 \q) = cosqG%- 2) (q) + S -^(P n _ 2 -l) / G^ 2) {q') sin 2q' dq 

7T 



G^ ( q ) = C0Sq G^>(q) + ^(/3„_! - 1) / G^' (q') sin 2q' dq 

n Jo 

I G [ N =0 (<l) = sing 
As in the case of bimodal disorder rewrite Eqs. (|33[) in a matrix form 

G(g, s) = F(g) + cos g M s G(g, s) -\ -BM a / G(g', s) sin2g' dq' 

t Jo 



,(n-l) 



sing 



(33) 



(34) 



where 



/ A) - i o 







B = 






02-1 






0„-2 - 









1 

0n-l - 1 / 



M. = 



/O 
1 
1 





Vo o 



1 \ 





1 
s 0/ 



and 



G(q,s) = 



G n -i(q, s) 
\G n (q,s) J 



F(?) = 



/ \ 




Introducing (as in Eq . (l29l) ) 



Q(s) = / G(<7, s) sm2qdq 
Jo 



we get for ([34 



i r 

Q(s) / dq sin q sin 2q 

t Jo 



I — cos q M s 



BM s Q(s) = dq sin2g 



I — cos q M s 



F(g). 



where by Z we have denoted the N x N unit matrix. 
It can be checked that BM S = M S B , where 



/ 0i - 1 
2 - 1 







03-1 



\ 








V o o o 



0n-l - 1 

n - 1 / 



0n = 00- 



Introducing the modified vector R(s) = B' Q(s) we rewrite Eq. (|3ip in the following form: 



- -l 1 
B / dq sin q sin 2q 



T S R(» 
with the right hand side 



The matrix T s in Eq. (|40[) can be explicitly written as 



I — cos q M s 



R(.s) = C(s) 



C(s) = / dq sin 2g 



I — cos q M s 



1 



a 



01-1 

a n-l GO 
Ora-2 



a 2 

Oi 



Ol 



1 



02-1 

a 2 



"2 
"1 



03-1 

"4 
"3 



F(g). 

Ctn-2 
a n-3 

1 



a 



0„-i - 1 

Ctn-l 



Ofi-2 



a 



0n-l / 
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and has disorder on its main diagonal only. The nonrandom elements a m are given by the following integrals 



a m (s) = — 



1 f n s sin q sin 2q cos" 



1 — s cos" q 



dq. 



(43) 



For n even the symmetry of the integrand implies that a m are zero for m odd. The structure of the matrix T s is 
therefore similar to that of Eq. ([31]) : it does not mix vectors having all but odd/even zero elements. In the following 
we call these subspaces odd and even sectors. It is important to note that C(s) belongs to the odd sector and 
therefore is affected by (3 m with m odd only. For n odd the integrand in the Eq. (|43|) has no symmetric properties and 
all elements a m are nonzero. In this case odd and even sectors are mixed by the matrix T s and the physical behavior 
of the polymer depends on all the values of the disorder potential. The behavior of the coefficients a m is depicted in 
FigE] separately for even and odd total lengths, n. In what follows we restrict ourselves to the case of odd values of 
n. This avoids the discussion of the somewhat pathological situation where the pinning transition is insensitive to a 
macroscopic number of values of the disorder. This situation arises due to the impossibility of return to the origin 
after an odd number of steps of our random walk and can be eliminated by a different choice of the step probabilities. 



Even n (n=200) 



Odd n (n=201) 




(a) 




FIG. 2: Behavior of the coefficients a m for even (a) and odd (b) values of total length n. 

The pinning transition point in the periodic potential /?' n ' = {Pi, ...,/?„} with the period n (n 3> 1), where n is 
assumed to be odd, is determined by the equation 



detf s {0 (n) } 



s=s tr = l 







(44) 



Let us stress once again that we are in the situation where the sequence {/3} consists of I = N / n copies of random 
subsequences {(3^} of length n (n 3> 1) each. In Appendix IBl we show that the pinning transition point in the chain 
consisting of £ (£ — > oo) copies of subsequences {f3^} {n — > oo) is the same as in the single subsequence (i.e. 
for I = 1) in the limit n — > oo. 



III. PERTURB ATIVE CALCULATION OF THE PHASE BOUNDARY 



using standard second 



In this Section we analyze the spectrum of the matrix T = T s given by Eq. 

s=s tT = l 

order perturbation theory. The disorder is supposed to be weak, i.e. the fluctuations of the diagonal elements in 
Eq. (|44|) are small compared to their mean value. 
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A. Non-random substrate 



Wc begin with the non-random situation and show that in this case Eqs. (|42j) -(|44 |) reproduce the results derived in 
Section [II Al Thus, we shall consider this non-random case as a reference state and develop a perturbation expansion 
with respect to this unperturbed state. In the absence of any disorder Eq. (|4"4")l reduces to 



detT s (/3) = dct 



( a - (0 - l)- 1 01 

On-i ao — (f3 - ly 1 



On-l 
dn-2 



01 



"2 



ao-(P-l)- 1 J 



where the non-random matrix T s {(3) is a special case of a Tocplitz matrix known as a circulant matrix [171 ] . We 
diagonalize it in the standard way: 



where A is the diagonal matrix of the eigenvalues 



-27rimk/n 



; m, k = 0, 1, Ti—l 



A,„ = a k e~ 2 ™ nk / n - (f3 - ly 1 m = 0, 1, ...,n - 1. 



Using the definition (|43| for the coefficients au and by virtue of (|46|) we have 



k=0 



1 f v ssmq sin2g cos 



n—l — k 



1 — s cos™ q 



dq 



e ~2iTimk/n _ Ijj _ 



At the critical value s = s tr = 1 the last equation becomes 



X m = 2e- 4 " m/ " - 1 - 2e- 4 ™' / " v / l-e 4 ™'/"- (/3 - l)" 1 . 
The transition point (3 tI is determined by the condition Xo(Pti-) = 0, yielding the criterion 

1 - (Ar - I)"' = 0; Ar = 2 

which reproduces the solution of the non-random problem considered at length in Section [II Al 



(45) 



(46) 



(47) 

(48) 
(49) 



B. Random substrate 



Suppose now that the energies u m = ln/3 m for m = 1, n are independent random variables distributed according 
to the law 



W + SCT n 



+ 1 with the probability ^ 
— 1 with the probability \ 



(50) 



and find eigenvalues and eigenvectors of the perturbed matrix T in the second order of series expansion in the strength 
of the disorder. Expanding the diagonal elements of the matrix T (Eq. ([42|) ) in e <C 1, we get 

Q, = a - (e"* - I)" 1 = «o - (e™ +,T * e - l)" 1 ; (i = 1, 2, n) (51) 
Expanding (|5ip , we arrive finally at the following expression for the diagonal elements Qi of the matrix T 

Q t = a - (e w - 1) - po + pia l (52) 
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where 



< 



le w (e w + l) o 
Po = --Hr. tt^-s 2 



2 (e w - l) 3 



(53) 



Pi 



(e w - l) 2 

The matrix T s {(3^ n '} up to the second order perturbation on s reads 



f a {pW} = A + B = f s (e w ) - p a I + 



I piCTi 
pi (72 



\ 



\ 





(54) 



where the matrices A = T s (e w ) — po! and B are correspondingly the unperturbed and perturbed parts of the matrix 
T (by / we have denoted the unit n x n matrix). 

Standard 2nd order perturbation theory leads to the following expression for the perturbed eigenvalue X' m : 

2 



(55) 

where v is the degeneracy of the eigenvalues and the eigenvectors v m are the columns of the matrix (I45[) . In what 
follows we arc interested in the computation of the largest eigenvalue Aq. 

For any particular sequence {(3^} the transition point is determined by the condition 



= \' {f3^} = 1 
Let us denote 



~ — 7 ~~P° 
e"- — 1 n 



In — 1 \ 9 / -i n-ln-1 



\k=0 



n* \ '■ — ' An — A,,, 

\ mVO k=0 k=0 



cr fc cr fc /e 



-27zm f (k— k' ) I n 



(56) 



n-l 
fe=0 



Ti — ± 71 — -L 

m'#0 A ° Am fe=0 fc=0 



(57) 



-2x771' {k — k') jn 



and expand the coefficients po and pi for small w in the vicinity of In 2 up to the 2nd order in w and e: 

75 

po = 3e 2 + 13e 2 (ln2 - 10) + y£ 2 (ln2 - w) 2 
Pi = 2e + 6e(ln 2 - io) + 13e(ln 2 - w) 2 
In the 2nd order perturbation series we keep only the terms up to the second order in w = In 2 — if and e: 



An 



-2w - 3w 2 - 3f 2 + (2s + 6ew)— + Ae 2 -^ = 



(58) 



(59) 



Solving the quadratic equation w(e) = (we have chosen the branch, on which w — > In 2 as e — > 0) and re-expanding 
the final expression in e up to the 2nd order, we get 



W 



^1 + le 2 (V3+4^ 



(60) 



Let us stress that the expression ([pTJ)) is valid for any distribution V{a} of random Ising-type variable <7j. The general 
form of the transition curve in the 2nd order approximation reads 



w 



Is 2 f 3^+4^1 



(61) 
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The terms proportional to (Si) vanish due to the symmetry er^ <-> — <7j. The calculation of the second order correction 
is straightforward though rather tedious and is given in Appendix [5J The final result for the averaged transition point 
w in the limit of large n reads 

* = y (62) 
so that the phase diagram is determined by the equation 

e 2 

w = \a2- — (63) 

This result of the perturbation theory for the phase boundary w(e) is shown in Fig[3]by dashed line. For e < 0.5 
there is a good agreement of the perturbative approach with the numerical data and results of our rcnormalization 
group computations of the next Section. 



IV. SIMPLE RG CONSIDERATION 



The location of the localization (wetting) transition can be found following a simple rcnormalization argument. We 
start with the master equation ([3]) in the Fourier space 

G N+1 (q) = cos qG N (q) + {fi N - l)-^ / dq' sin 2q' G N (q') (64) 



and iterate it twice 

G N+2 (q) = cos 2 qG N (q) (65) 

+ (fa - l) SingC ° Sg f dq' sm2q'G N (q') + (p N+1 - l)^ f dq' sin 2q' cos q' G N (q') (66) 
^ Jo 11 Jo 

+ (0 N -l)(0 N+l -l)—*. dq' sin 2<z' sing' / dq" G N (q"). (67) 
n Jo JO 

The last term in the r.h.s. vanishes due to the orthogonality. The second and third term describe processes of two-step 
arrival to x = 1. We note that in the long-wavelength limit q — > they arc identical as seen by replacing the cosines 
by unity. This corresponds to neglecting variation in Gn{x) on scale of the lattice spacing and will be justified later. 
Combining these terms we obtain an expression similar to the small q version of Eq. (|64| with modified random part 

G N+ 2(q)-G N (q) = -q 2 G N (q) + (f3 N + f] N+1 -2)^- dq' sin 2q'G N (q') (68) 

n Jo 

and corresponding to the double time interval. By this procedure the random potential has been normalized to the 
arithmetic mean of two subsequent terms: 

Pn/2 = 2 ( ' 

The probability distribution of the renormalized disorder j3' is related to that of f3 as follows. Let P{0) be a proba- 
bility distribution of (3j and P'((3') the one of the new variables. The corresponding Fourier transforms (characteristic 
functions) £(A) and £'(A') are related as 

e'(A') = e 2 (A'/2). (70) 
The fixed point of this transformation is the exponential function 

£'(A)=£(A)=exp(iA</?» (71) 



which is just a consequence of the Central Limit Theorem. 
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Therefore the disorder (3j = e Uj can be just replaced by their mean value ((3j) = (e Uj ) in the long-wavelength limit. 
The crucial observation is that (e Uj ) e~'"', so fluctuations of Uj modify the localization criterion. For disorder 
discussed in the previous section (see (|50|)). The mean value of (3j is given by 

fa) =e w coshe (72) 

so it depends on both parameters, the mean value, w, and fluctuations, e. The localization criterion reads from (|17[) 

j3 = e w coshe = 2 =>■ io = In 2 - lncoshe (73) 



The lowest order in expansion of this expression in powers of e yields the result of the perturbation theory ()63|) . The 
graphical representation of the RG phase boundary (|T3"|) is shown in FigJ3]by solid line. 
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0.0 



0.5 



1.0 
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2.5 
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FIG. 3: Phase diagram for quenched disorder in the plane (e,w). The renormalization group computations and perturbation 
theory are depicted by solid and dashed lines correspondingly. The results of numerical averaging are shown for odd n 
(n = 1001). 



Finally, we observe that from the expressions (j9|)- ([T0|) it follows that the bound state contribution to the Fourier 
transform 



G N (q) 



(3-2 sing _ N 
2(3 l-s cosg S ° ' 



so 



(3 



(74) 



is strongly peaked around q = close to the transition point (3 = 2, This justifies replacing cosg by unity in the 
integrals in Eq. (|65| . 



NUMERICAL ANALYSIS OF THE DEPINNING TRANSITION 



To visualize the eigenvalues of the matrix T , consider the case of absence of any disorder, i.e. take in Eq. ((42|) 
the homogeneous sequence {(3^} with (3j = (3 for all j = l,...,n. In Fig 0] we have plotted the eigenvalues A m 
(m = 1, ...,n) of the matrix T in the complex plane for two values w = In 2.1 and w = In 2.0 (transition point) for 
even (a) and odd (b) values of n. For even values of n the spectrum of the pure matrix T is doubly degenerate, while 
for odd n this degeneracy disappears. In what follows we consider odd n only. 

We now consider disorder generated by the distribution (|50p and find numerically the location of the pinning 
transition. The general procedure is very simple: we fix some value e, then generate ensemble of random sequences 
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FIG. 4: The eigenvalues A m of the matrix T in the complex plane (ReA,ImA) for: (a) even n (n = 200) and (b) odd n {n = 201) 
for two different values of w: w = In 2.0 (transition point) and w = In 2.1. 



p( n ) f r om the distribution (|50[) and find for each sequence such a w at which the real eigenvalue of the matrix T 
crosses zero implying the condition 



lim detf /?„} = 



(75) 



The obtained critical value is then averaged over realizations {/3^} of quenched disorder. We take the size of the 
n x n matrix T equal to n = 1001. In Fig[3]we show the points corresponding to the averaged phase boundary for 
the ensemble of quenched sequences in the space of the parameters e,w. In the next section we discuss the 

obtained numerical phase diagram. 



VI. CONCLUSIONS 



In this paper we have discussed the problem of wetting of a one-dimensional random walk on a half-line x > in 
a short-ranged periodic potential located at the origin x = 0. When the period of the potential becomes equal to the 
length of the random walk trajectory, then one could conjecture that the trajectory under consideration is interacting 
with the quenched random potential. In other words, we approach the random potential, increasing more and more 
the period of the substrate potential. 

The most attention is paid to the question concerning the location of the transition point averaged over all equally 
weighted realizations of the disorder. We have shown that for any specific disorder realization, the transition point 
can be obtained from the condition that the determinant of the matrix (|4^|) (with the entries (|4"3"|) ) is equal to zero 
at s = 1 — see Eq. (|4"4"|) . The period of the substrate is given the size n of the matrix (jH?)) which remains to be very 
large but finite in our consideration, and the fugacity s is conjugate to the number of copies, £ = N/n, of rj-pcriodic 
potential (so, N is the total length of the substrate). Setting the fugacity s in ([4^1) - (|^|) to its limiting value Str = 1 
amounts to the consideration of I — > oo copies of the system. 

One may worry that in the thermodynamic limit the transition point in the periodic sequence of I (£ — » oo) 
subchains each of length n (n» 1) differs from the transition point in a single (£ = 1) subchain (of the same primary 
structure) of length n (n — > oo). However this is not true and in Appendix [Bl we show that the transition point is not 
sensitive to the sequence of thermodynamic limits and for n — > oo is independent on number of copies, £. 

We have compared the results of our numerical simulations performed in a wide range of disorder strengths with the 
results of the perturbation theory for weak disorder, and the results of a simple renormalization group consideration. 
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The results obtained allow us to conclude that for weak disorder e < 0.5 all approaches (including the perturbation 
theory) arc in a good agreement. For sufficiently high disorder strengths e > e pcr t ^0.5 the second order perturbation 
theory fails, while the results of the renormalization group always agree with the numerical simulations — see the FigEJ 
This point deserves some special discussion. Namely, recall that results of the RG computations claim the existence 
of the averaged Boltzmann weight which governs the behavior of the system under renormalization. This fact signals 
the existence of an effective preaveraged (with respect to the disorder) annealed system with averaged Boltzmann 
weight (e Um ). This statement is consistent with the conclusion of works @, but contradicts the statement that the 
disorder is marginally relevant formulated for the first time in Q. 

Let us note at the very end that the advantage of our matrix approach is its transparency. We have made no any 
uncontrolled assumption and have reduced very complicated initial problem with many degrees of freedom to the 
well posed problem of finding eigenvalues of some random matrix with a relatively simple structure. The localization 
criteria detT{/3i, ...,/3/v} = for a random sequence {(3\, ...,/3/v} in the limit N — > oo is exact for any strengths and 
any distributions of (3i (i — 1, N). 

The results obtained in our work allow us to state that even for strong disorder the shift of the averaged pinning 
point of the random walk in the ensemble of random realizations of substrate disorder is indistinguishable from the 
pinning point of the system with preaveraged (i.e. annealed) Boltzmann weight. 
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APPENDIX A: PERTURBATION THEORY FOR MATRIX T s 

The terms in Eq. (foTj) are 

n-l 

n 



I 1 

= -£<**) 

n n — ' 

fe=o 

i(^ 2 ) = iEE(^) (ad 



k=a k'=o 

k=0 k=0 

The computation of (Si) and (Sf) where averaging (...) is performed over the symmetric distribution V{<j\, is 
straightforward: 

-(Si)=0; ^{S!) = - (A2) 

Using (|48|) . we obtain for the unperturbed part A of the matrix T (Eq. (|54|) ) the following expression for the 
eigenvalues 

A m = 2e- imm/n - 1 - 2e- 4 ™/Vl - e 4 ™"/» - (e w - l)" 1 - Po (A3) 

Taking into account that 



y J l _ e 4W« = .Asin^e 1 ^™/"- 1 / 4 ), 
V n 

we can substitute (|A3[) into last line of (|A1|) . Thus we get 

J_ (5 2 ) = Iy 1 =iV (2 - 2e - 4 ™"'/" + 2e--( 1 / 44 - 3m V»)./ 2 sin— ^ (A4) 

Since the eigenvalues are distributed symmetrically with respect to the real axis (see FigJ^i), the sum ^2 
takes only the real values. After some algebra, we arrive at the final equation for the expectation (S2)' 



m An— A„ 



1 9 («-i)/2 1 -cos^HL + A /2sin^^cosf^ + %) 

\{S 2 )= 2 - V " V L= ( A5 ) 



m '#o sin 2^ f 1 - 2,/2 sin sin (^™1 + + 2 sin 
In the limit n — > 00 the sum in (|A5[) can be replaced by the integral 



1 /" 1/2 1 - cos Attx + V2sin27ra; cos (37ra; + 

- <^> « 2 / 7 ^ d - dx = 1 (A6) 

Jo sin 27rx ( 1 - 2 V2 sin 2ttx sin (3nx + f ) + 2 sin 27 



Substituting (|A2|) and (| A6|) into [) we arrive in the limit n — > 00 at the desired equation (p2"|) . 
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APPENDIX B: INDEPENDENCE OF PHASE BOUNDARY ON NUMBER OF PERIODS I IN THE 

THERMODYNAMIC LIMIT 

It follows from the general definition of the generating functions, that in order to consider the phase transition in 
an individual period of the substrate potential (|32[) . i.e. for I = 1, we should shift the fugacity s (conjugated to £) 
away from its marginal value s = Str = 1 (corresponding to i — ► oo) to some value sg lying in the interval < sg < s tr . 

The number of periods, £, is controlled by the fugacity si. So, to compare the transition points in chains of different 
periods, we should normalize the Boltzmann weights of chain of I periods per a weight of a a single period, si (for 
I = 1). Dividing the Boltzmann weights in the matrix T s (see Eq . (|4"2"|) i by si, we get: (j3 — 1) — » (/?* — l)/s^- Now we 
could investigate how the normalized transition point, /3j* r , depends on the typical fugacity sg. It is easy to see that in 
the absence of any disorder (i.e. for e = 0) the eigenvalue Aq (see Eq. (jT7|) for m = 0) is given by the following integral 



The equation Xo(si) = in the limit n — > oo determines the position of the normalised transition point, P^ r (si). 
One can easily verify that in the limit n — > oo the integral (|B1[) is independent of Sf and, hence, /3j r (sf) = 2 for any 
< si < Str- The same conclusion holds for e ^ in the matrix T s : the transition point in the limit n — > oo for any 
random primary sequence is independent of the effective fugacity sg corresponding to finite number of periods, I. 

So, we conclude that in the thermodynamic limit the transition point in the periodic sequence of i (£ — > oo) random 
subchains each of length n (n 1) coincides with the transition point in a single (£ = 1) random subchain (of the 
same primary structure) of length n (n — > oo). 

This conclusion has very transparent physical sense. Namely, consider from the very beginning the random walk 
with periodic boundary conditions in the space. The simplest way is to suppose that the first (t = 0) and last (t = N) 
steps of the random trajectory are always attached to the random substrate at the point x = 0. For such a periodic 
system the location of the transition point is exactly given by the equation 



Hence, the transition point is not sensitive to whether the terminal step of the random walk is attached to the surface 
or not. This is evident in the localized regime where the fluctuations of mean-square end-to-end distance are constant 
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(Bl) 



lim detT{/3 1; /3 2 ,..,^v} =0 



(see Eq.©). 



